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Abstract: In this paper, the investigation and detailed modeling of a cascaded Raman laser, 
operating in the midwave infrared region, is described. The device is based on 
silicon-on-insulator optical waveguides and a coupled resonant microcavity. Theoretical 
results are compared with recent experiments, demonstrating a very good agreement. 
Design criteria are derived for cascaded Raman lasers working as continuous wave light 
sources to simultaneously sense two types of gases, namely C2H6 and CO2, at a moderate 
power level of 130 mW. 
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1. Introduction 

For over two decades compact, broadly tunable, energy efficient midwave infrared (MWIR) and 
longwave infrared (LWIR) sources and devices have been the topic of active research [1]. Historically, 
the need for sources operating especially in the 3-5 j^m and 8-12 jam atmospheric transmission 
windows has been primarily driven by military applications such as wind light detection and ranging 
(LIDAR), and IR countermeasures (IRCM). However, in recent years such sources have also found 
use in a wide array of applications ranging from purely scientific uses, such as ring down and Fourier 
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transform infrared (FTIR) spectroscopy, to clinical and industrial uses such as tissue ablation and 
hydrocarbon detection [2]. In addition, the growing interest for industrial uses such as hydrocarbon 
detection from vehicle, oil fields, and industrial smoke stacks has recently induced the research to 
increase its efforts to optimise and study lasers for mid infrared gas sensing. 

Laser-based gas sensing is attractive because it can provide a way to achieve highly sensitive, 
real-time, in situ detection of various gases. The key to its success is the development of compact, 
reliable, and cost-effective laser light sources operating in the infrared region, where various gases 
exhibit a large number of strong absorption lines originating from fundamental, rotational and 
vibrational absorption [3,4]. Currently, near-IR distributed feedback (DFB) laser diodes developed for 
telecommunications are used for gas sensing because they satisfy all these requirements. However, 
there is a strong need to develop mid-IR laser sources because the absorption intensities of most gases 
are well larger in the mid-IR than near-IR by a factor 100-10,000. As shown in [3], a large number of 
very important gases for industrial, environmental and safety needs, such as HC1, CH 4 , C2H6, CO2, 
NH 3 , N 2 0, S0 2 , H 2 0, can be detected using lasers operating in the mid-IR wavelength range 2-10 jam. 
Thus, the use of mid-IR lasers is expected to greatly increase the sensitivity of gas sensing and reduce 
the optical path length and system sizes. 

Recently, mid-IR laser sources based on difference frequency generation in quasi-phase-matched 
(QPM) LiNbC>3 have been studied as promising candidates for these applications because they can 
provide continuous wave (CW) mid-IR light in the wavelength range 2-5 jam at room temperature. 
However, in spite of their excellent sensitivity, the low conversion efficiency of the conventional 
QPM-LiNbC>3 devices has limited their use because large, expansive, high power lasers must be also 
used to achieve a reasonable amount of mid-IR optical output [3]. Moreover, different structures of 3.5 |um 
small wavelength distributed feedback quantum cascade lasers have actually been developing to 
achieve low threshold mid-IR sources for gas sensing [5]. A recent work proposed in [6] has shown the 
possibility to use an interband cascade laser as light source to detect gases in the range "k = 3.6 + 4.3 jam. 

Another approach recently proposed to realise mid-IR light sources is based on Stimulated Raman 
Scattering (SRS) effect. In fact, one of the major advantages of Raman lasers is their ability to generate 
coherent light in wavelength regions that are not easily accessible with other types of lasers [7]. To this 
aim, silicon is a particularly suitable material for Raman lasers operating in the near and mid-IR 
regions because it can guarantee a very good trade-off between low cost and high performance. As 
well described in [8,9], silicon represents the ideal platform for Integrated Optics and Optoelectronics 
because the quality of commercial silicon wafers driven by Microelectronics industry continues to 
improve while the cost continues to decrease. Moreover, the compatibility with silicon integrated 
circuits manufacturing and silicon Micro-Electro-Mechanical Systems (MEMS) technology is another 
important reason for this interest in Silicon Photonics. As a transmission medium, silicon has much 
higher nonlinear effects than the commonly used silicon dioxide, in particular the Raman effect. In fact, 
Raman gain has been successfully exploited in fiber amplifiers and lasers, but usually several 
kilometres of fibre are required to create a useful device. Fortunately, the gain coefficient for SRS is 
approximately four orders of magnitude larger in silicon than in silica. Additionally, 
silicon-on-insulator (SOI) waveguides can confine the optical field to an area that is approximately 100 times 
smaller than the modal area in a standard single-mode optical fiber. These two circumstances allow the 
fabrication of efficient Raman-based centimeter-scale integrated optical devices to be achieved. 
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Starting from 2002, several experimental and theoretical studies based on this effect have been 
proposed in literature, such as Raman amplification in SOI waveguides [10-15], Stokes and 
anti-Stokes Raman conversion [16-19], cross phase modulation-based interferometer switch [20] 
two-photon absorption [21,22], lossless modulation [23] to point-out the real efficiency to use SRS 
effect on the Silicon platform. More recently, efficient Raman lasing in silicon in near-IR (NIR) region 
has been experimentally and theoretically demonstrated [24-28], showing great potential for realising 
low-cost, compact, room-temperature lasers in MWIR [1,29-33]. 

Therefore, in this work we theoretically analyze the possibility to realise a Raman cascaded laser 
for simultaneously sensing two different gases having their absorption peaks in the wavelength 
range 3-5 jam. Our choice is motivated by two concurrent aspects. First, the cascade Raman emission 
can result attractive in the range 3-5 jam, since two photon absorption (TP A) and free carrier 
absorption (FCA) effects are completely eliminated. Second, the cascade Raman laser could represent 
an efficient answer to the drawbacks of interband diode lasers and quantum cascade lasers. In fact, the 
ability to induce the cascaded lasing by means of SRS effect can potentially produce several 
wavelengths in the range 2-5 jam, at the same time too long for interband diode lasers to be reached 
due to Auger recombination, and usually too short for quantum cascaded lasers owing to the 
finite-conduction-band offset at room temperature [34,35]. However, two works have also 
demonstrated the possibility to realise light sources by means of interband diode lasers and quantum 
cascaded lasers emitting around 3.8 jam [6,36]. 

This paper is organized as follows. In Section 2 we derive the mathematical model to study the 
nonlinear effects in a resonant microcavity coupled to an external waveguide, following a different 
approach from literature [30]. The proposed modeling includes all nonlinear effects involved in the 
integrated structure without any a-priori assumption, including SRS, Self-Phase-Modulation (SPM) 
and Cross-Phase-Modulation (XPM) effects as induced by Kerr nonlinearity, nonlinear birefringence 
effect as induced by the modal birefringence and walk-off The model considers the interaction 
between pump pulse and first-order Stokes wave, as well as excitation of higher-order Stokes waves, 
mismatch between the input beams and microcavity resonance wavelengths, and coupling mechanism 
between microcavity and bus waveguide. In Section 3 a number of numerical results are presented, 
including comparisons between our theory and some experiments in literature on CW microcavity-based 
cascaded Raman lasers for validating our approach. Moreover, a systematic study of SOI waveguides 
modal properties, the polarization dependence of coupling factors of both pump and Stokes waves and 
the laser performance are investigated at MWIR for the first time, to the best of our knowledge. Finally, 
Section 4 summarizes the conclusions. 

2. Device Modeling 

In this section, a very accurate physical model proposed in our previous works [37-39] is 
generalized to analyse the Raman lasing effect in a SOI microcavity resonator working at mid-IR. The 
model is based on a set of partial differential equations for nonlinear coupling among pump, first-order 
and higher-order Stokes waves inside the microcavity, considering both polarization states. 

Modeling of cascaded Raman laser in mid-IR region has been already proposed in literature [30]. 
However, differently from that approach where the source is obtained by placing the silicon waveguide 
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between two highly reflecting elements such as Bragg gratings or dielectric coatings, we consider the 
racetrack resonator coupled externally with the bus waveguide as laser architecture. This fundamental 
difference induces to extend the model in [30] to include in the physical analysis the resonance 
condition mismatch of input pump and Stokes waves as originated inside the resonator, as well as the 
effect of photon decay time induced by the coupling process. In addition, differently from other 
experimental and theoretical works presented in literature [30-32], we also consider in this paper as the 
waveguide size influences the modal birefringence, the group velocity mismatch and, thus, the laser 
performance. 

In our analysis we assume the architecture as sketched in Figure 1, where the input pump (S p ) is 
injected in the resonant microcavity by evanescent coupling from an external SOI bus waveguide. As it 
was already experimentally demonstrated [26,27,31], the most efficient architecture for the laser cavity 
is represented by a racetrack resonator. In fact, it leads us to control the proper combination of both 
pump and Stokes coupling factors by changing the coupling length L coup , gap G and polarization state. 

Figure 1. Schematic architecture of a cascaded racetrack-resonator SOI Raman laser. 




The scheme of Figure 1 also clarifies the concept of a cascaded racetrack-resonator SOI Raman 
laser. The pump coupled into the laser cavity generates optical gain by SRS inside the silicon 
waveguide at the first-order Stokes wavelength, which is 15.6 THz red-shifted from the pump. The 
gain increases with the pump power, and lasing threshold is reached when the optical gain equals the 
total cavity loss. The optical power of the first-order Stokes wave starts to increase inside the cavity, 
inducing a pump depletion and eventually gain for the second-order Stokes signal. This cascaded 
process can continue and generate higher-order Stokes lasing at longer wavelengths. It is worth to 
outline that the cascaded racetrack-resonator SOI Raman laser shown in Figure 1 represents the laser 
head for the gas detection equipment. In fact, an experimental setup to demonstrate the gas detection 
should involve driver, laser head, filters, gas cells, photodetectors and electronic system for signal 
processing [3]. Thus, the input pump laser modulated by the driver produces the mid-IR outputs as 
described above, that are collimated by means of filters into the gas cells and detected by 
photodetectors. Finally, the photocurrents are processed by the electronic circuit. 

Hereinafter, we assume a SOI waveguide as sketched in Figure 2, having rib total height H, slab 
height H s and rib width W. Without any lack of generality, the electric field inside the microcavity is 
predominantly assumed as a single transverse mode. This condition can be satisfied by choosing 
appropriate rib sizes to meet the single-mode condition. 
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Figure 2. Cross-section of SOI waveguide. 

y 




SiQ 2 



In a single-mode SOI waveguide two propagating modes are typically confined, one quasi-TE 
(dominant x-component of electric field) and one quasi-TM (dominant y-component). Our coupled-mode 
approach describes the power transfer among pump wave (p), first-order Stokes wave {s\) and higher-order 
Stokes waves. Two higher-order Stokes waves (52,53) and both polarizations are considered. 

Under assumption of translational invariance along the propagation direction (z), due to large radius 
R of racetrack resonator (see Figure 1), the electric field in the single mode SOI waveguide can be 
written using the variable separation principle as E(x,y,z,t) = C-F{x,y)A{z,t)e ikz , where k is the 

propagation constant and F(x,y) is the optical mode distribution in the waveguide cross section x - y 

(solutions of Helmholtz wave equation), C is a normalization constant, c=\j j \F (x,y)\ 2 dxdy > an d 

A(z,t) is the slowly-varying wave amplitude. Thus, without any leak of generality, the total electric 
field inside the SOI waveguide can be written as: 

I q»A „(*, t)F, m {x, y)e^—>) + X ^(z, t)F Xn (x, y)e ^-^ 



E(x,y,z,t) = x 



+y 



m n 



(1) 



where the meaning of superscripts and subscripts is as follows: 1 = p {TE \ 2 =p (TM) , 3 = s^ 1 "' , 4 = s l 



(TM) 



5 = s 2 (TE) , 6 = s 2 (TM) , 7 = s 3 (TE) , 8 = s 3 (TM) . In this equation, summations in x and y are relevant to 
quasi-TE and quasi-TM modes, respectively, subscripts m , n , r and q designate all possible cavity 

longitudinal modes for pump and first, second and third Stokes waves, respectively, co p>m (co s] , n , co s2r , a> S 3, q ) 
is the resonant angular pulsation of the pump (Stokes) mode inside the cavity and k t m (i = 1,2) is the 

pump propagation constant under resonance condition. 

As detailed in [38], it is clear that only a few terms are to be considered in the summations of 
Equation (1) under external excitation. In fact, if the input pump (S^)is launched in the bus waveguide 
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with angular frequencies a> p , the first-order Stokes wave is originated at co s =co p - Q.r, being Q. R = 15.6 THz 
the Raman frequency shift in silicon. In this condition, only those modes whose resonant angular 
frequencies are closer to co p and co s i will give a contribution to Equation (1). Similar considerations 
hold for higher-order Stokes waves. These modes are characterized by an longitudinal order, namely 
in , n , r , and q , evaluated by means of the resonance condition as m~ co p n eff p L cavity> j(2nc) , 

n~co sX n eff A L cavity ,j{lnc) , r ~ co sl n effsl L cavity f{2nc) and q ~ co s3 n^ s3 L cavity /(2xc) , being 

n e ff,p{n e ff, s \, n e ff, S 2, «e#>3) the effective index of pump (Stokes) wave inside the racetrack resonator, c the 
light velocity in vacuum, and a> p , co s i, co S 2, a> S 3 the angular frequencies of pump, first-order, second 
order and third-order Stokes waves, respectively. 

Moreover, the previous consideration holds if the cavity free spectral range (FSR) is larger than the 
input pulse bandwidth, Aa> pu i se . For example, assuming a Gaussian input pulse with full-wave-half- 
maximum (FWHM) time width T F whm, the condition FSR > Aa> p i U s e (pulse bandwidth) gives 
L cavit y < ( 2 xcT FWHM )/(2yj2ln2n efLp ) , being L cavity the cavity perimeter [38]. Then, by assuming 

Tfwhm =100 ps, the cavity length required to consider only two resonant modes for the pump wave 
and aligned as quasi-TE and quasi-TM modes can be estimated to be smaller than 2.8 cm, that is a 
reasonable condition in some integrated structures. Now, by following the procedure outlined in [37,38] 
and assuming the nonlinear contributions to V NL as a small perturbation of refractive index, we have 
obtained the following coupled equations for pump waves: 

8A 8A l 8 2 A l 8 3 A 11 2 

V ^ + lf + 7 2 V ^V"6^ ^ 

+7' 2v g,^ |4| 2 4 +i' 2v g ,^3 |4| 2 4 +J 2v g ,£^ A |4f 4 + 2 7\A^,A4*4^ <A,J+ * 4 ~ V * 3)Z 

+ 2jv g ^ AD ,AAAe J ^ 

1 CO , |2 1 CO 1 ,2 7], 

-^^£R—ftAH A £-^^ R —feAH 4+7 — v sA 

si si cavity 

£ = 1,2 (2) 

being £=1,2 for quasi-TE and quasi-TM polarization, respectively. 

Then, we have for quasi-TE and quasi-TM first (k = 3,4) and second order Stokes {k= 5,6) pulses: 
8A„ OA . 1 „ d 2 A„ 1 „ d 3 A^ .,_ x . 1 1 



v 



dz + ^ + j v g ,^p2, K ^-\, K ^^ = j{co K -co K )A K --— A^+jlv^M 4 

I 4+> g ,Jvl4| 4+7 2 w-,v|4| 4 + 2^^/^141 4 
+v g , K \g R ,J^i |4 + if 4 + / 2 v,^^, /I+1 ,,4 + i^4^ ( "" +,+ ^ F """ JZ + / 2 v„^^ +1 ^4^ + i4^' ( "" + ^"" +1 " Fjz 

- 2 jr j 1 I |2 1 II 2 

+ i v gAwA4'i e " * ~~ 2~ v S'* r ^'* : / |: '/' +4 1 4+4 1 4 ~~ 2~ v S'' i '^' /i '' c '/'^+ 5 14+5 1 4 



, - ? l A A* A p j( k " +k ^- k *- k ^ z 4. ,-o v i- A A* A J^ +k i^- k '- k ^ z 



/r = 3,4,5,6 (3) 



and for quasi-TE and quasi-TM third order Stokes (p = 7, 8) waves: 
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v 



p = 7,S (4) 



8z dt 2 g ' p ' zp df 6 g ' p ' >' p df Jy p 2t 

h 

I |2 I i2 II 2 II 2 

+ /v 2r , L4, ^4 + jv y AAA A + iv y \A \ A + jlv y AAA A 

J g,P 'P,5\ 5\ p J g,P'P,6\ 6 1 p J g,p'p,p\ p\ p J g,P' P,i \ S P 



+ \ v g , P S R J P , 5 14| 2 A p +-v g „g R jf p6 \A 6 f A p + 72v gi / pA6i ^X4 e ' ( " 6+% ^"* 5)Z 

V g,p p,6,5,f 5 6 f J g,p p,p,C,C C P C 

In Equation (2) subscript v = 2,1 holds when E,= 1,2, respectively. In Equation (3), if *r= 3,4, 
^ = ^i ' ®* = c ° s \w ^ = 1 an d ^ = 4,3, while if k = 5,6 , = <y j2 , co K =co slT ju = 3 and v = 6,5. 

Moreover, K - g R — - . In Equation (4), it holds £ = 8,7 when q = 7,8 respectively, and 

g R3 -g R — — . In conclusion, Equation (2) takes into account the time-space evolution of both 

quasi-TE and quasi-TM pump pulses. Similar evolutions are described by Equation (3) for first and 
second order Stokes pulses of both polarizations, and by Equation (4) for third order Stokes pulses of 
both polarizations. In equation system (2)-(4), the asterisk denotes the complex conjugate. The partial 
differential equation system (2)-(4) clearly represents a generalization for mid-IR laser applications of 
the model proposed in [37,38]. The terms (co pm -co p ), (co, lw -co sl ), (co s2T -co s2 ) and (co s3Ij -co s3 ) 

designate the mismatch from the resonance condition of input pump and Stokes waves originated 
inside the resonator, respectively. The term r^with E, = 1,2 (r„r p ) represents the overall photon decay 
time of the pump (Stokes) wave inside the cavity. It can be expressed as a function of the following 
two terms: 

— = — + (5) 

where Xi^, and x c ^ represent the decay time related to total linear losses inside the cavity and coupling 
between the optical resonator and external bus waveguide, respectively. As detailed in our previous 
work [37], these terms can be calculated as: 

J = J_A 

being v g ^ the pump wave group velocity, a!°* s , the overall linear loss coefficient related to 
propagation, bending and scattering loss contribution, and 77? the coupling factor, defined as the power 

fraction of the pump wave leaving the resonator due to the directional coupler. Similar definition holds 
for r r and t p . 

In Equations (2)-(4), g R = 3cQ s ijUoX R (Q R )/(n e ff iP n e f£ S \) [2] is the Raman gain with //„ the magnetic 

susceptibility, n e ff iP , and n^r s i the effective refractive index of the pump and first-order Stokes waves, 
respectively, and z R (Qr) the Raman-resonant susceptibility, defined as in [16]. Moreover, +g R 
(positive sign) designates the SRS effect, while -g R (negative sign) determines the depletion effect due 
to the energy exchange induced by SRS for n -th order Stokes wave, if (n + l)-th order Stokes wave is 
simultaneously present in the relevant equation. 



A 1 cavity / r\ 

T l'=l \- Z" c = — r- (6) 
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Moreover, the coefficients ytj = nioofjc and = ri2(ofijlc take into account SPM and XPM effects 
as induced by Kerr nonlinearity, being c the light velocity and ni the nonlinear refractive index [39,40]. 
Terms kijxi = nicofij^i/c, i,j = 1,...,8, depend on birefringence effects. The overlap integrals fij and 
fijxi are given by: 



jF^y)] |F,(x,,y)| dxdy ^ g 

Jjl^ (/x^JJ|f, (x,j)| dxdy 

+00 

\ J F * F *j F k F i dxd y 

fi,j,k,i = ~f / N/ \/ \ Uj,k,l = 1,...,8 

J J |i^| 2 dx£/y J jr ;'| J J \ F kfd x dy J J \F^dxdy 

v —GO / \ — 00 / V — 00 / V — 00 y 

In particular, = A eff . represents the effective core area of the optical mode relevant to the i-th 

pulse, i = 1,...,8. Terms kijjj are usually dominant {kixjj > hj.ki f° r ' ' ^ j> k ^ 0 an d represent the 
coherent coupling between the polarization components for waves at the same frequency, giving 
degenerate four- wave mixing [41]. Rigorously, additional terms should be introduced in Equation (4) 
to take into account XPM and FWM effects between third-order Stokes waves and pump waves. 

However, since these terms can be considered as negligible they are not included in Equation (4). In 
fact, we can assume FWM effect between third-order Stokes waves and pump waves as negligible due 
to a too large phase mismatch, while the XPM effect between the third-order Stokes waves and the 
pump waves can be considered negligible with respect to XPM effect between third-order and 
second-order Stokes waves, due to the strong pump depletion effect induced by the first-order Stokes 
wave. In addition, in continuous-wave regime all terms related to SPM and XPM effects do not 
influence strongly the laser behavior since Raman susceptibility dominates over Kerr susceptibility in 
SOI technology. 

For completeness, group velocity dispersion (GVD) and third-order dispersion (TOD) effects are 
also included in the model, as indicated by the terms proportional to fii.i and /? 3>i . In any case, these 
effects are removed from physical description in case of CW Raman laser applications. On the 
contrary, with pulsed Raman lasers GVD and TOD effects are not negligible for very short pulses. 
However, in this case, input pump pulse Tfwhm is much shorter than z c ^, then the pump energy inside 
the resonator cannot be enhanced. Thus, as demonstrated in [38], pulsed excitation in resonator-based 
Raman lasers should require T F whm» r c ,^to induce the SRS effect above the threshold. 

Finally, it is worth to note that TPA and FCA effects are not included in our equations due to 
excellent transmission of silicon in mid IR because of the absence of TPA for wavelengths longer 
than 2.2 um [2]. Thus, cascade Raman emission can result very attractive in the range 3-5 um, since 
TPA and FCA effects are completely absent. 
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3. Numerical Results and Discussion 

3.1. Comparisons with Experimental Measurements 

To test the analytical formulas and physical assumptions for Raman effect into optical microcavities 
under CW operation, we have compared our numerical results with some experiments proposed in 
literature. A very interesting set of comparisons with experimental results involves the CW cascaded 
Raman laser based on SOI resonator. The architecture used in the experimental setup [31] is the same 
as Figure 1, where the microcavity is constituted by a SOI racetrack resonator working with a pump 
beam at 1.55 um. The silicon waveguide cavity used in the experiment was fabricated on a SOI wafer 
with a 1 um buried oxide layer, using CMOS compatible fabrication processes. The proposed device 
has rib width W= 1.5 um, height H= 1.55 um, and etch depth 0.7 um. The coupler gap is 0.7 um, and 
the coupler length ranges between 900 and 1,100 um to obtain the desired coupling coefficients for 
pump and signal wavelengths. The total length of racetrack cavity is 3 cm and the bend radius is 400 um. 
The bus waveguide is 1.6 cm long and it is connected to the ring cavity by means of the directional 
coupler. Finally, the measured optical parameters include ai oss = 0.2 ± 0.05 dB/cm, r e ff = 0.4 ns, 
gR= 10.5 cm/GW, and fi TPA = 0.5 cm/GW, being ff PA the two-photon absorption coefficient. In order 
to achieve cascaded lasing, it is crucial to design and fabricate the directional coupler to have the 
proper combination of coupling factors for pump and Stokes waves. To this aim, we have assumed the 
same coupler geometry as in [31] and the relevant coupling factors have been evaluated by means of 
a 3D CAD tool based on beam propagation method (BPM) [42]. Thus, our simulations are not 
conditioned by any arbitrary choice or fitting parameters. 

It is very important to consider also TPA and FCA effects in Equations (2)-(4), since the laser 
pump is tuned at 1.55 um. Thus, according to our previous work [37], we include the contribution 
a\ FCA) , i = 1,2,... 8, to the total losses, as induced by free carrier density change generated mainly by 

TPA of pump pulse. In addition, a\ FCA) is evaluated according to Soref s relationship [43] as: 



a, (FC4) =8.5-10- 18 



( 



X. 



v 



1.55 



f 



AiV, +6.0-10 



-IS 



1.55 



2. 



>2 



1.55 



N c (7) 



where N c = AN e = AN/, is the density of electron-hole pairs generated by TPA process. The coefficient 
oq = 1.45 x 10~ 17 cm -2 [16] is the FCA cross section measured at X = 1.55 um, and A; is the relevant 
mode wavelength (either pump or Stokes wave). Finally, the rate equation governing the free carrier 
dynamics into the waveguide core is required together with system (2)-(4), given by [12]: 

at r eff 2na> p V 1 / 

where z e ff is the relevant effective recombination lifetime for free carriers, h is the reduced Planck 
constant, and ff PA is the coefficient of carrier generation induced by TPA process of the pump beam. 
In turn, this nonlinear absorption effect requires the inclusion of contributions 0 .5v ' g ,£P TPA j \^Ai 
VgJ TPA fKA A i A * and v g ,ffi TPA fpJAi 2 A p in Equations (2), (3) and (4), respectively. 

Figure 3 shows the laser output power (Stokes power at the bus waveguide end) versus pump input 
power. The black and red markers represent the experimental data for first and second order Stokes 
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waves, respectively. In experiment, the pump wavelength is at 1,550 nm and the first and second lasing 
wavelengths are measured at 1,686 nm and 1,848 nm, respectively. 

Figure 3. Comparison between experimental data in literature and modeling of this paper 
in terms of CW cascaded Raman emission versus input pump power (X p = 1.55 urn, 
X sl = 1.686 urn, X s2 = 1.848 (am). 




Pump power (mW) 



The solid lines designate our numerical results as evaluated by solving the coupled equations 
proposed in the model. Our numerical calculations demonstrate the same group velocity inside the 
resonator for pump and Stokes waves, and thus the partial differential equations (2)-(4) can be 
transformed in ordinary differential equations by simply writing [38]: 

dA v , dA B s dA B , 

^l = ^i + v — 5± (9) 
dt 8t 8 dz 

where v gtP = v gsS \ = v gj a = v g = 8.2645 x 10 7 m/s. The plots show a very good agreement with 
experimental data in terms of threshold values, output powers above threshold, external efficiencies 
and output saturation. In particular, the good agreement with first Stokes output saturation confirms 
how the Raman cascaded lasing is well described by the mathematical model proposed in this paper. 

3.2. Waveguide Optical Properties in Mid-IR 

To the best of our knowledge, systematic investigations have been not yet presented in literature for 
optical properties of SOI waveguides in mid-IR region, although this is a topic of increasing interest [44]. 
The goal of this sub-section is to determine the MWIR optical characteristics of SOI waveguides in 
terms of birefringence and group velocity mismatch. 

As it is well known, a single-mode waveguide can in general support two modes polarized in 
orthogonal directions. Under ideal conditions, a mode excited with its dominant polarization, i.e., in x 
direction (quasi-TE mode), would not coupled to the mode with orthogonal dominant y-polarization 
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state (quasi-TM). However, in real waveguides random variations of cross section shape and 
stress-induced anisotropy result in a mixing of two polarization states. Thus, the two modes exchange 
their powers in a periodic fashion as they propagate inside the waveguides with period 
L B = Xlun™ • It is worth to note that the modal birefringence could induce the nonlinear 

birefringence [39], as considered in the model through coefficients kij,kj- In fact, since the electric field 
associated with an arbitrarily polarized optical wave can be written as in Equation (1), the nonlinear 
part of induced polarization can be given in similar form but each component depending on third-order 
susceptibility (Kerr effect) [41]. In fact, since the third-order nonlinearity involves in general 
interaction among four optical waves, following the algebra details in [41] each component of the 
induced polarization nonlinear part can be given as the sum of a number of terms proportional to «2 
and to the products between four electric field components. Some of these terms are responsible for 
SPM and XPM effects, while the remaining terms in ktj&i are related to FWM. In addition, as 
explained in [41], significant FWM effect occurs only if the phase mismatch nearly vanishes. This 
should require matching of both frequencies and wave vectors. In case of CW Raman laser applications, 
the nonlinear birefringence could become evident as a rotation of the polarization ellipse [41]. This effect 
can be neglected by appropriately choosing the cavity length. 

Figure 4 shows the modal birefringence ( An eff = n™ - n"i ) spectra for different values of r = HJH. 

In particular, in our simulations we have assumed a rib total height H = 2.2 um, and a rib width W 
of 1.5 um (solid lines) and 2 um (dashed lines). Hereinafter, calculations of waveguide optical 
properties in mid-IR region have been carried out by full-vectorial finite element method (FEM) [45]. 

Figure 4. Modal birefringence spectra for different values of r and W. 
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The plots show that for each value of the rib width, the modal birefringence changes sign with 
increasing r from 0.2 to 0.5, increasing with wavelength. In addition, W= 2 um leads a reduced value 
of \An e fj\ for each value of r and in all wavelength ranges to be achieved. The curves of Figure 4 are 
useful to calculate the beat length Lb and thus the minimum cavity length needful to avoid the 
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nonlinear birefringence effect. In fact, it occurs L cavity »LB for long waveguides with large 
birefringence, then the terms containing kijxi often change sign and the total phase contribution due to 
birefringence averages out to zero [39]. In contrast, this contribution is not negligible if L cavi y < L B (kajj), 
as it occurs in short waveguides with moderate birefringence. Using the simulations proposed in 
Figure 4, we can observe that for W = 1.5 um, r = 0.2 and operating wavelength 5.5 urn, the modal 
birefringence has maximum value, inducing a beat length of Lb = 55 um. Thus, a cavity length much 
larger than 55 um should be required to average to zero the contributes depending on kij,kj- In 
particular, it is very interesting the case with W = 2 um and r = 0.3, where the curve shape is flat 
around zero in 2.5-3.5 um wavelength range. For example, for a wavelength of 3.5 um, we obtain 
Lb -6.92 mm. Therefore, it is evident that waveguides with small birefringence require large cavity 
lengths in contrast with the general requirement of small chip area occupation. However, the cavity 
length has to be determined to simultaneously satisfy the condition L cavit y » Lb and the condition able 
to induce SRS effect into the cavity. 

Another important optical characterization of SOI waveguides in mid-IR region can be made in 
terms of group velocity mismatch between quasi-TE and quasi-TM modes. Figure 5 shows the walk-off 
parameter defined as d TMTE = v~ 1 (TM)-v~ l (TE) versus operating wavelength, for different values of 

W and r. Similarly to simulations in Figure 4, the rib total height H has been assumed of 2.2 um. 
Figure 5. Walk-off parameter spectra for different values of W and r. 
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As a general trend, the spectral curves of Figure 5 show that the walk-off absolute value increases 
with increasing the wavelength. However, three cases are worthy to be pointed out. The first 
corresponds to the waveguide width W= 1.5 um and r = 0.5, where an absolute maximum takes place 
around the operating wavelength of 5.22 um. The second case is related to a SOI waveguide with 
W = 1.5 um and r = 0.3, where a zero crossing occurs at a wavelength of 4.82 um. In this case, it is 
evident how two optical pulses, polarized in orthogonal states and injected into the waveguide, will 
propagate with the same velocity. Finally, it seems to be very interesting the case for W = 2.0 um and 
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r = 0.2. For this waveguide, we can observe a flat shape around zero level in the wavelength 
range 2.5 ^ 3.5 \im. It is important to outline that this last case induces both weak modal birefringence 
and very small group velocity mismatch between quasi-TE and quasi-TM modes. 

3.3. Design Guidelines for Cascaded Raman Lasers 

Recently, efficient Raman lasing in silicon has been demonstrated in NIR region [26,27], showing 
great potential for realizing low-cost, compact, room-temperature lasers even in MWIR region [31]. 
Such lasers are highly desirable for many applications ranging from trace-gas sensing, environmental 
monitoring and biomedical analysis, to industrial processes control and free-space communications. To 
the best of our knowledge, only a few papers are reported in literature on MWIR Raman lasers design. 
In fact, in the most of experimental and theoretical papers, as in [30,31], any comment about influence 
of waveguide size and mode polarization on the laser performance are not reported. Thus, the goal of 
this sub-section is to suggest a number of numerical simulations as main guidelines to design both 
laser optical waveguide and directional coupler, with the aim to optimize the lasing effect. 

Without any loss of generality, hereinafter we focus our attention on the design of a Raman laser 
source for ethane gas sensing. In Figure 6 the ethane spectrum is sketched, characterized by an 
absorption peak centered at 3.3485 \im. Thus, we assume the first-order Stokes wavelength of 
l s \ = 3.3485 fxm for ethane gas detection. This implies a laser pump tuned at X p = 2.8519 (am 
(antimony-based laser), while the higher-order Stokes waves are originated through cascaded Raman 
effect at/l S 2 = 4.0544 (j,m and/U = 5.1376 fxm, respectively. 

Figure 6. Ethane spectrum around 3.35 \im. 
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In the following analysis, we guess that the pump wave is aligned with quasi-TE polarized mode (as 
usually occurs in most experimental set-ups). Under this assumption, Equation (2) holds only for £, = 1. 
In addition, since quasi-TM polarization is filtered out, the slowly-varying wave amplitude Ai is zero 
and then all terms in ktjxi and proportional to have to be set to zero in Equation (2). However, 
Equations (3)-(4) hold their validity because a quasi-TE pump pulse could in principle generate Stokes 
waves of both polarizations. 
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Figure 7. Walk-off parameter between first-order Stokes and pump waves for different 
values of// and r (X p = 2.8519 um, A s i = 3.3485 um). 
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For the following discussion, we assume that the first-order Stokes wave grows-up being mainly 
aligned as a quasi-TM mode. This assumption will be verified in the next sub-section, due to smaller 
optical area and thus larger Raman modal gain of quasi-TM over quasi-TE mode. Under this 
assumption, we can evaluate as the walk-off parameter between quasi-TM first-order Stokes wave and 
quasi-TE pump is influenced by the waveguide sizes. Then, Figure 7 shows this walk-off parameter 
versus waveguide width, for different values of rib total height H and r. 

As a general trend, Figure 7 shows the walk-off increase with decreasing H for any value of both r 
and waveguide width. Moreover, for each value of H the curves show a monotonically decreasing 
shape as a wavelength function in case of r = 0.5, while assuming the opposite trend for r = 0.3. In 
case of r = 0.4, the curve is monotonically decreasing or increasing for H = 2.6 um and H = 1.8 um, 
respectively. In addition, a very interesting case is represented for r = 0.3 and H = 2.6 um, where a 
zero crossing at W = 1.5334 um occurs. Thus, a waveguide designed with this optimal size can 
guarantee the rigorous absence of walk-off effect from a physical point of view, according to relation (8). 

As demonstrated in [38], the walk-off effect could represent the limiting factor for pulsed Raman 
lasers. In fact, the walk-off length L w = 7o/|c?tm,te| could be shorter than the cavity length 
(T 0 = 7WWL665 for Gaussian pulses). Then, SRS is limited by the group-velocity mismatch and 
occurs only over distances z~L w , even if the cavity length L cavity , » L w . At the same time, the nonlinear 
effects such as SPM and XPM become important because of the relatively large peak powers, 
considerably affecting the evolution of both pump and Raman waves. Thus, our optimal waveguide 
designed in this sub-section definitely avoids the walk-off problem in any case, including pulsed 
operation. Then, SRS effect induced in the resonant microcavity could be only limited by the 
enhancement factor of the pump wave. 
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Figure 8. Modal birefringence for pump and first-order Stokes waves versus waveguide 
width for different values of H (X p = 2.85 19 um, l s \ = 3.3485 um). 
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Figure 8 shows the modal birefringence for pump and first-order Stokes waves versus waveguide 
width, for different values of H and r = 0.3. The plot shows zero crossings for both pump and Stokes 
waves only for H < 2.2 um. In particular, the birefringence free waveguide width decreases with 
decreasing H. The curves shown in Figure 8 can be used to estimate the beat length Lb- In particular, 
for the optimal waveguide (r = 0.3, H = 2.6 um and W = 1.5334 um), we calculate L B (pump) = 214.42 um 
and L B (\ st Stokes) = 174.4 urn. 

Thus, an optical resonator with L cav ity » 214.42 um has to be used to minimize the nonlinear 
birefringence effect. In this condition, it could be assumed that two Stokes waves of the same order but 
aligned in orthogonal polarizations are simultaneously amplified in a SOI waveguide by the same 
pump pulse, without any reciprocal interference. Their different behaviour will only depend on 
effective area and coupling factor. 

Further, we show in Figure 9 both GVD and TOD coefficients as a function of waveguide width for 
both pump and Stokes waves. In the simulations we have assumed r = 0.3, H= 2.6 um and W= 1.5334 um, 
and both polarizations. It is interesting to note how GVD coefficients assume a negative value for both 
pump and first-order Stokes waves and both polarizations. In addition, GVD and TOD coefficients are 
larger than those obtained in SOI waveguides operating in near-IR region [39], opening important 
challenges for nonlinear signal processing in mid-IR. 
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Figure 9. GVD and TOD coefficients versus waveguide width for pump and Stokes waves 
{k p = 2.8519 urn, X sl = 3.3485 urn). 
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Now, both optical waveguides and minimum cavity length are designed and estimated. Another 
very important aspect is to find the design guidelines for the directional coupling. In fact, as shown by 
a formula derived in [38], the output efficiency and the lasing threshold for first-order Stokes depends 
strongly on the coupling factors related to the waves propagating inside the optical resonator. In 
particular, to reduce the threshold level for the first-order Stokes wave, it is essential to maximise the 
enhancement effect into the racetrack resonator. By defining the enhancement factor as T = \A^ max /S^, 
being A^ max the pump amplitude maximum inside the cavity and Sg the input pump peak outside the 
cavity for a fixed polarization state, the design of architecture in Figure 1 requires the condition T » 1 . 

Due to the absence of detrimental effects induced by FCA and TPA contributions in mid-IR, it is 
possible to estimate a priori, with good approximation, the condition to maximise the enhancement 
factor. In fact, by considering a simple linear coupling mechanism and neglecting all nonlinear effects, 
it is possible to demonstrate that Y factor depends essentially on the ratio K = r c ^ between total 
loss and pump coupling decay times. Following the standard conventions, under-coupling is denoted 
by K < 1, over-coupling by K > 1 and critical coupling by K = 1. This last condition means vanishing 
bus waveguide transmission, and thus it represents the optimal condition to maximise the pump energy 
inside the optical resonator, minimizing the Raman lasing threshold. 

The previous guideline is rigorous if we are mainly focusing on first-order Stokes wave emission. In 
the case of cascaded Raman laser [31], it is generally appropriate to design the coupler close to the 
critical coupling for the pump wavelength; close to zero coupling for the first Stokes wavelength 
(achieving high intracavity power that generates Raman gain at the second-order Stokes wavelength) 
and low coupling for the second-order Stokes, in order to obtain at the same time both low lasing 
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threshold and enough output power extraction from the resonator [31]. In our work, the design criteria 
for directional coupler are substantially different. In fact, with the aim to design a Raman laser for 
sensing of two different gases, polarization selectivity and output power extraction of the same order 
of magnitude for the first and second-order Stokes waves are required. 

A number of 3D simulations based on BPM show a good trade-off between previous requirements 
when we consider a gap G = 0.8 urn and a coupling length L coup = 1,486 um, being the coupling 
factors 3.13% and 2.47% for first and second-order quasi-TM Stokes waves, respectively. In turn, 
these values guarantee low threshold and similar output power extraction for both Stokes waves. 
Figure 10 shows the power exchange in the directional coupler versus coupling length for the pump 
wave tuned at 2.8519 um and aligned as a quasi-TE mode, for G = 0.8 um. 

The plot also shows a quasi-TE pump wave coupling factor of 5% at L coup = 1,486 um. Thus, 
assuming the total propagation loss inside the cavity a'° 0 ' ss = 5.0657 m -1 (0.22 dB/cm) [39], the critical 

coupling condition for the pump (77? = L cavity ,a'°' ss ) requires a coupling factor ranging from 5.06% 

to 20.26% as the cavity length increases from 1 to 4 cm. Therefore, hereinafter L cavity = 1 cm is 
assumed. In Table 1 the coupling factors for Stokes waves in directional coupler are summarized. 

Using the waveguide sizes and coupling factors, we can progress to the cascaded Raman laser 
simulation in the range 2.85 4.0544 um. According with some experimental works (i.e., [2,29,31]), 
we have considered in our numerical investigations silicon chip orientation [10 0] and Raman gain of 
gR = 9 cm/GW. However, to the best of our knowledge, experimental indications are not mentioned 
about any Raman gain anisotropy with respect to pump beam polarization. In addition, we assume 
n 2 = 5 x 10 " 5 cm 2 /GW, J3 TPA - 0 and negligible FCA effect [2]. 

Figure 10. Power exchange in the directional coupler for pump wave (X p = 2.85 19 um). 
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The pump is assumed aligned as a quasi-TE mode with n e jf y \ = 3.2911, A e ff t \ = 2.4585 um 2 and 
v gi i = 8.3424 x 10 7 m/s. The optical waveguide parameters, related to Stokes waves and obtained by 
means of full-vectorial FEM method [45], are summarized in Table 2. 
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Table 1. Coupling factors for Stokes waves. 
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Table 2. Optical waveguide parameters for Raman laser simulations. 
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To complete the cascaded Raman laser design, some comments about the absorption losses in the 



mid-IR region are needed. The total propagation loss inside the cavity a^°'J is strongly influenced by 

the material absorption losses that can become the dominant contribution in the mid-IR region. In fact, 
SOI technology in the mid-IR is usually limited by the silica that present an absorption loss coefficient 
less than ldB/cm over 0.25-3.6 (am, except for a spike of 10 dB/cm in the range 2.6-2.9 (am [1]. 

Thus, according with results proposed in [1], the total absorption loss coefficient (cc^' } ) for SOI rib 

waveguide can be estimated as: 

«r=rXr+r Si02 <f (io) 

being a { ^" lk) and oc^ lk) the bulk absorption of silicon and silica, respectively. The terms Y Si and 
T SiQ ^ are guided-power fractions into silicon and silica layer, respectively. Equation (10) suggests that 

it is possible to guarantee a modal low optical loss coefficient if SOI rib waveguide is designed with H s 
large enough to reduce the field tails inside the silicon oxide layer. As mentioned before, the critical 
wavelengths for the cascaded Raman laser as proposed in this work could be X p = 2.8519 jam, and 
X S 2 =4.0544 jam, where silicon oxide presents large absorption values cc { s b ^ k) ~10dB/cm, 

and ~7dB/cm, respectively. However, the silicon layer presents a very low absorption value, i.e., 
a ibuik) _ 0.001 dB/cm for both wavelengths, and only a very small fraction of optical power is confined 

inside the silicon oxide. Thus, our estimation of using Equation (10) and summarized in Table 3 

demonstrate that the optimistic goal of a*° 0 ' ss = 0.22 dB/cm can be really achieved with the SOI 

waveguide designed and proposed in this work. 
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Table 3. Total absorption coefficients. 
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Figure 1 1 shows the output of the cascaded Raman laser as designed and operating in continuous-wave 
regime. In our analysis, the laser output is calculated at the bus waveguide end as: 

ll.-iiJ'S'" ' (11) 

where P s is the power inside the cavity for the specific Stokes wave and L Bu s represents the total bus 
length (see Figure 1). 

In the plot we can observe the output power for first and second-order Stokes waves considering 
both polarizations. Several comments are worthy to give. First, the quasi-TE second-order Stokes 
wave presents a large value of lasing threshold, thus it cannot be lasing with an input power less 
than 150 mW. On the contrary, the quasi-TM second-order Stokes presents similar lasing threshold 
and external output efficiency, independently from the polarization state of exciting first order Stokes. 
It seems very interesting the first-order Stokes wave lasing. In fact, in this case the coupling factors for 
quasi-TE and quasi-TM polarizations do not represent a discriminating factor for the threshold level, 
since it has the same order of magnitude. In any case, above threshold (equal for both polarizations) 
the quasi-TM mode presents an output power level and an external efficiency larger than that for 
quasi-TE mode, as induced by smaller core area and larger coupling factor. Thus, the quasi-TM 
polarization is well suitable for the simultaneous lasing of the first two Stokes waves, allowing 
detection of two gases. 

Moreover, Figure 11 demonstrates an output power of 7.136 mW for both first and second-order 
Stokes waves for an input power of 131.25 mW. Under this condition, the time dynamics for the 
designed cascaded Raman laser is sketched in Figure 12. It is worth to note as the time reference 
moves with the waves, as indicated by relationship (9), and the designed cascaded Raman laser 
operates in continuous-wave regime. The previous simulations demonstrate that not only the cascaded 
Raman lasing can be achieved in SOI waveguides working at mid-IR, but also that its performance are 
well suitable for gas sensing environmental applications in terms of reduced power consumption, 
stability, wavelength tunability and polarization selectivity. In addition, our theoretical investigations 
have demonstrated the possibility to obtain output powers larger than 7 mW, much better than in other 
approaches, i.e., 0.2 mW as mentioned in [3], with consequent positive impact on relevant electronic 
signal processing and device sensitivity. 
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Figure 11. Laser output versus pump power for different combinations of polarization states. 
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In particular, the Raman laser is designed to detect two different types of gases, i.e., C2H6 and CO2 [3]. 
It is evident that, in spite of an increasing input pump power, the third-order Stokes wave lasing at 
X S 3 = 5.1376 um could be also obtained, giving for example the possibility to detect a third gas, such 
as NO [3]. 

Figure 12. Time dynamics for both pump and Stokes waves (X p = 2.8519 um, 
l,i = 3.3485 um, X s2 = 4.0544 um). 
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4. Conclusions 

The generalized model presented in this paper allows to accurately predict the time dynamics of 
pump and Stokes waves in cascaded Raman lasers based on a SOI microcavity resonator and operating 
in mid-IR region. In addition, the theoretical study has allowed to demonstrate high performance of 
cascaded Raman laser in terms of reduced power consumption, stability, wavelength tunability and 
polarizations selectivity, offering unique advantages as mid-IR optical sources and thus competing 
with complex and bulk solid-state laser systems. Finally, a detailed design of a cascaded Raman laser 
for gas sensing of C2H6 and CO2 has been developed under CW operation. This approach could be also 
used for gas sensing at longer mid-IR wavelengths with improved sensitivity. 
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